Magnetostriction in the Bose-Einstein Condensate quantum magnet NiCl2-4SC(NH2)2 
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The quantum magnet NiCl2-4SC(NH2)2 is a candidate for observing Bose-Einstein Condensation 
of spin degrees of freedom in applied magnetic fields. An XY antiferromagnetic ordered state occurs 
in a dome-shaped region of the temperature-field phase diagram between Hci = 2.1 T and Hc2 
= 12.6 T and below 1.2 K. BEC corresponds to the held-induced quantum phase transition into 
the ordered state. We investigate magnetostriction in single crystals of this compound at dilution 
refrigerator temperatures in magnetic fields up to 18 T, and as a function of magnetic field angle. 
We show that significant changes in the lattice parameters are induced by magnetic fields, and argue 
that these result from antiferromagnetic couplings between the Ni spins along the tetragonal c-axis. 
The magnetic phase diagram as a function of temperature, field, and field angle can be extracted 
from these data. We discuss the implications of these results to Bose-Einstein Condensation in this 
system. 
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I. INTRODUCTION 

In recent years, there has been a surge of interest in the 
topic of quantum magnetism in condensed-matter sys- 
tems. Nature and human ingenuity have provided us 
with a rich variety of spin structures including reduced- 
dimensional ladders, chains or planes, dimers, frustrated 
spins, and single-molecule magnets. These materials all 
share the common trait that quantum effects, such as 
spin fluctuations and quantized spin levels, play a signif- 
icant role in shaping the ground state and the physical 
properties of the system. The resultant behavior can be 
complex and challenge our understanding, as well as pro- 
vide potentially useful applications. 

In particular, there have been a number of attempts 
recently to observe the quintessential quantum ground 
state Bose-Einstein Condensation (BEC) in quantum 
magnets. BEC was first observed in dilute gases of *^Rb 
atomSfi leading to a nobcl prize being awarded in 2001. 
It turns out that a form of BEC can also be observed 
in quantum magnets, e.g. crystalline lattices containing 
spins. BEC in quantum magnets was first predicted to 
occur in 1991 by Ian Affleck^. In the past few years, sev- 
eral reports of BEC in real spin systems have been pub- 
lished including TlCuClg,-^'^ BaCu2Si06,-^>^ CsCuCU,-^-^ 
and NiCl2-4SC(NH2)2-^'^ 

The bosons that condense in these systems are not the 
atoms, but rather Cu or Ni spin degrees of freedom (also 
referred to as triplons or magnons in the literature.) The 
BEC corresponds to the onset of magnetic order induced 
by magnetic fields. Thus, the tuning parameter for induc- 
ing condensation in spin systems is not the temperature, 
but the magnetic field. 



In this paper, we focus on the organic quantum magnet 
NiCl2-4SC(NH2)2 (DTN). It is the first compound based 
on Ni (spin S = 1) instead of Cu (S = 1/2) and the first 
organic compound in which BEC has been observed.— 
We present an extended introduction to the topic of 
Bose-Einstein condensation in DTN and then discuss new 
magnetostriction data for this compound at low temper- 
atures. 

The energy level diagram of the Ni system differs from 
those of the previously studied Cu systems. The Ni 5 = 1 
spin triplet is split by single-ion anisotropy into a, Sz — 
ground state and 5"^ = ±1 excited states with an en- 
ergy gap of £» ~ 10 K.iaJJ.a2 The 5"^ = ±1 levels are 
dispersed by the antiferromagnetic coupling J, e.g. the 
spins can have different energies depending on their ori- 
entation with respect to each other. Thus these levels 
become broad bands, as shown in Fig. [2] where the bot- 
tom of the band is the AFM k = (tt, tt, tt) vector. When a 
magnetic field is applied along the tetragonal c-axis, the 
Zeeman effect lowers the Sz = 1 level until the bottom 
of the band becomes degenerate with the Sz = ground 
state at Hd (see Fig. [2]). Between Hd and Hc2 the 
mean-field ground state becomes a linear combination of 
all three energy levels until finally at Hc2 the gap reopens 
now with the Sz = I level as the ground statci^ 

At zero field, the gap between the Sz ~ and the 
S'z = ±1 levels precludes any magnetic order. However, 
between H^i and Hc2, XY antiferromagnetic order oc- 
curs at low temperatures. The region of field-induced 
antiferromagnetic order between Hd and Hc2 and below 
the maximum = 1.2 K is shown in the phase diagram 
in Fig. [TJ At Hd, the spins order antiferromagnetically 
in the a-b plane, perpendicular to the applied field direc- 
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FIG. 1: Temperature T - Magnetic field H piiase diagrarr 
for H\\c determined from specific lieat and magnetocalorie 
data.i^ The magnetization vs field at 16 mK is overlayed onto 
the phase diagram." 
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FIG. 2: Energy level diagram of DTN. The = ±1 states 
are shown schematically as broad bands due to antiferromag- 
netic dispersion. The spins (blue arrows) form a disordered 
spin liquid with Sz = below Hci, then order in a canted 
AFM structure between Hci and Hc2- Finally, above Hc2 
the spins polarize along the field direction. The ground state 
wave function > in these three regions is also indicated. 



tion. As the field increases, the spins cant along the field 
direction until finally at Hc2, the spins polarize and the 
magnetization saturates. The magnetization along the c 
axis is shown overlaid on the phase diagram in Fig. [1] 
We should note that this description only holds for fields 
applied along the tetragonal c-axis. For fields in the a-b 
plane, the Sz ~ groundstate mixes with a linear com- 
bination of the Sz = ±1 excited states and the resulting 
new set of eigenstates move apart with field, thus there 
is no level crossing and the system does not order. 

Inelastic neutron scattering measurements in zero 
field^ indicate that the antiferromagnetic coupling is 
strongest along the Ni-Cl-Cl-Ni chains along the tetrag- 
onal c-axis, with the AFM exchange parameter Jc ~ 



FIG. 3: Unit ceU of tetragonal NiCl2-4SC(NH2)2 showing Ni 
(red) and CI (blue) atoms. Remaining atoms have been omit- 
ted for clarity. 



1.74(3) K. The crystal structure of DTN is shown in Fig. 
[31 The coupling along the a-axis (and cquivalently the 
b-axis), Ja = Jb = 0.17(1) K, is significantly weaker, and 
no couplings were found to within the experimental res- 
olution along the (1,1,1) direction. Thus the magnetic 
structure of this compound consists of semi-ID chains of 
Ni atoms that order three-dimensionally below 1.2 K in 
fields. 

In DTN, the Bose statistics come about not because 
the spins are S = 1. In fact, theoretically BEC could be 
observed in a purely S = 1/2 magnetic system. Rather, 
the Bose statistics result from the fact that only the spins 
play a role in the Hamiltonian, and the atoms themselves 
are confined to lattice sites. The spin operators on dif- 
ferent sites commute with each other, leading to Bose 
statistics. 

The spin S ~ I system can be mapped into a gas 
of semi-hard-core boson o^^i^^'^^ with no more than two 
bosons per site. The \Sz = I >, \Sz = > and 
\Sz = —1 > states are mapped into the states with 
two, one and zero bosons respectively. The transverse 
component of the Heiscnberg interaction is mapped into 
a hopping term while the longitudinal or Ising compo- 
nent becomes a nearest-neighbor repulsion. The single- 
ion anisotropy D generates an on-site repulsion and the 
magnetic field is mapped into a chemical potential. The 
resulting Hamiltonian in the particle representation is an 
extended Hubbard model for the semi- hard-core bosons. 
In the particle representation, the phase for H < Hd 
corresponds to the Mott insulator (one boson per site). 
When the chemical potential (magnetic field) reaches a 
critical value H = Hd , the density p of bosons (magneti- 
zation) starts to increase continuously and the additional 
bosons relative to p = 1 form a BEC state at low enough 
temperature. This state is the particle version of the XY- 
antiferromagnet. When the chemical potential reaches a 
second critical value H = i/c2, all the sites are occupied 
by two bosons (p = 2) and the system enters a second 
Mott insulating phase that is the particle version of the 
fully saturated spin state. 

The tetragonal crystal structure of DTN creates an 
approximate uniaxial symmetry (U(l) symmetry) of the 
spin system. This guarantees that the antiferromagnetic 
order for fields between Hd and Hc2 is XY-like, meaning 
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that there is no energetically favorable direction for the 
spins to point in the plane perpendicular to the applied 
field. The order parameter therefore acquires not only 
a magnitude but also a phase and the resulting Bose- 
Einstein Condensate exhibits spontaneous phase coher- 
ence. 

There are several caveats to the idea of BEC in quan- 
tum spin systems. First of all, the uniaxial spin sym- 
metry of the Hamiltonian is approximate. The square 
lattice of the crystal introduces anisotropy in the a- 
b plane as a result of spin-orbit coupling and dipole- 
dipole interactions.— However, these effects generally oc- 
cur at lower energy scales (often fiK) and can thus be 
neglected at the temperatures of tens to hundreds of mK 
at which these quantum magnets are studied. In par- 
ticular, no significant Dzyaloshinskii-Moriya interactions 
have been observed in electron-spin-resonance measure- 
ments of DTN^ 

Thus, the Bose-Einstein condensation picture in quan- 
tum magnets is valid for the temperatures at which these 
compounds are studied, and more importantly it provides 
a way of understanding the observed physical behavior. 
The thermal phase transition in this system belongs to 
the d = 3 universality class of an XY antiferromagnet, 
where d is the number of spatial dimensions. However, 
the field-induced quantum phase transition belongs to 
the D = 5 {d — 3, z = 2) universality class where z is 
the dynamical exponent. 

One of the key experimental signatures of the Bose- 
Einstein Condensation quantum critical point is the 
temperature-dependence of the critical field Hd near 
T = 0. For a 3-D BEC, HdiT) oc where a 
3/2. In contrast, the prediction for an Ising magnet is 
a = 2, and for a 2-D BEC, a = 1 ) (T corresponds 
to the Berezinskii-Kostcrlitz-Thouless transition temper- 
ature in this case). In DTN, Hci{T) as T ^ has been 
investigated in detail by specific heat and magnetocaloric 
effect measurements)^ These data show that the expected 
power-law behavior Hci{T)-Hci{G) oc T^/^ for 3-D BEC 
is approached as T ^ Oi^ 

In this work, we present magnetostriction measure- 
ments of DTN performed in a dilution refrigerator in a 
20 T magnet system. The original motivation for these 
experiments, as mentioned in a previous work;^ was to 
account for a discrepancy between the observed and the 
predicted phase diagram. A semi-classical spin-wave ap- 
proach predicted that Hc2 ~ 10.8 T, in contrast to the 
observed value of iJc2 = 12.6 T. In fact, it turns out that 
this discrepancy results from the presence of strong quan- 
tum fluctuations that limit the validity of a spin- wave ap- 
proach. Taking fluctuations into account, the calculated 
value of Hci is very close to the experimental one. Quan- 
tum Monte Carlo calculations are able to model both the 
T—H phase diagram and the magnetization as a function 
of field very wellii^ The magnetostriction is therefore not 
required to account for any deviations in the observed 
results. 

Nevertheless, upon measuring magnetostriction in this 



compound we have found significant field-induced lattice 
distortions in DTN, and these indicate that magnetoe- 
lastic coupling plays a role in this system. We find field- 
induced changes in the lattice parameters up to 0.025%, 
which is a relatively large effect. We also map the mag- 
netic phase diagram as a function of magnetic field angle 
from the tetragonal c-axis. 



II. RESULTS AND DISCUSSION 

Magnetostriction measurements of single crystals of 
DTN were performed in a titanium dilatometer.'' The 
sample was mounted on a titanium base plate using Si 
vacuum grease. Length changes were monitored capaci- 
tativcly with a Bc-Cu spring-mounted titanium tip and 
the sample was protected from the tip with a .003" Ti 
foil. The dilatometer was mounted in a plastic rotator 
inside the '^He-'*He mixture of a top-loading dilution re- 
frigerator, in a 20 T superconducting magnet system at 
the National High Magnetic Field Laboratory in Talla- 
hassee, FL. Temperature was monitored with a ruthe- 
nium oxide thermometer, and capacitance changes were 
measured using a digital capacitance bridge operating at 
5 kHz (Andcen-Hagerling model 2700A). Data is shown 
above _ff = 1 T, since flux jumps in the NbaSn supercon- 
ducting magnet affect the data at low flelds. 

New magnetization data for DTN is also presented. 
The longitudinal magnetization for H\\a was measured 
in a vibrating sample magnetometer (VSM) at 0.5 Kf^i 

Magnetostriction data at 25 mK is shown in Fig. [4] 
for magnetic flelds along the crystallographic c- and a- 
axes (top and bottom flgures, respectively). For both 
fleld orientations, magnetostriction was measured along 
the a- and c-axes of the crystal, e.g. both perpendicular 
and parallel to the applied fleld. 

For H\\c (Fig. S] top,) changes in the length of the 
sample along the c axis, ALc, clearly show features at 
the magnetic phase transitions at Hd and Hc2- ALc 
is nearly constant below Hd and above -ffc2, but varies 
dramatically in the region of antifcrromagnctic order be- 
tween Hd and Hc2- Between Hd = 2.1 T and 5 T, Lc 
shrinks, and between 5 T and Hc2 = 12.5 it expands, fl- 
nally saturating at Hc2 ■ This behavior is rather startling 
since the magnetization for H\\c increases roughly lin- 
early with field between Hd and 7?c2- The longitudinal 
magnetization for H\\c is shown for comparison in Fig. 
O We interpret this nonmonotonic behavior in ALc as 
arising from competing effects of antifcrromagnctic and 
ferromagnetic components of the canted order. At Hd 
the magnetic system orders purely antiferromagnetically 
in the a-b plane. Thus, by shrinking the distance between 
Ni atoms along the c-axis, the system can increase Jc and 
lower its magnetic energy. However, with increasing field, 
the spins cant along the field direction, increasing the fer- 
romagnetic component. At some point it becomes more 
energetically favorable for the system to expand the c- 
axis thereby lowering the antifcrromagnctic coupling Jc- 
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FIG. 4: Normalized percentage length change %AL/L as a 
function of magnetic field measured along the crystallographic 
c-axis (solid blue lines) and a-axis (dashed red lines). Data 
is taken at 25 mK. The magnetic field is applied along the 
c-axis (top), and a-axis (bottom). The inset to the top figure 
is zoomed in on the feature at Hci in %ALa/ La- A straight 
line has been subtracted from the data for clarity. 



Finally, above i?c2 the spins are saturated along the ap- 
plied field direction and remains roughly constant. 
The magnitude of the overall change in Lc between Hd 
and Hc2 is ^ .025%. Such a distortion would be consid- 
ered large in a metal and would most likely correspond to 
a structural phase transition. However, since DTN is a 
soft organic compound, a .025% distortion of the lattice 
parameters is reasonable. 

The a-axis magnetostriction AL^ for H\\c, shown as 
the red line in Fig. 0] (bottom), is smaller ALc and Hd 
is barely distinguishable at 2.1 T. The smaller magne- 
tostriction along a makes sense since Jo << Jc- What 
is interesting is that AL^ doesn't follow the same non- 
monotonic behavior as ALc- One might expect at first 
glance that ALa should echo the same behavior as ALc, 
just with a smaller magnitude. In reality, as the mag- 
netic field becomes larger, ALa doesn't reverse directions 
like ALc does as the spins cant, but rather continues to 
decrease until it saturates near Hc2- This odd behavior 
could be due to the Poisson effect, e.g. volume-conserving 
forces typically present in tetragonal crystals cause La to 
contract as Lc expands. Thus at high fields, the a-axis 
lattice parameter responds more strongly to the distor- 
tions in the c-axis parameter than to magnetic forces from 
the antiferromagnetic coupling Ja- 

For the other magnetic field direction, H\\a, no mag- 
netic order has been observed, and the system behaves 
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FIG. 5: Longitudinal magnetization of DTN for -ffjjc and 
H\\a. The H\\c data is taken at 16 mK in a force 
magnetometer ,- 
T = 0.5 K. 
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FIG. 6: Temperature T - Magnetic field H phase diagram 
for -ffjjc determined from magnetostriction, with the phase 
diagram from specific heat and magnetocaloric measurements 
shown for comparisoit^. The magnetostriction phase diagram 
(squares) is shifted to slightly lower fields, which may result 
from a slight pressure applied by the dilatometer. 



paramagnetically. The magnetostriction shown in Fig. [4] 
(top) reflects this, and evolves smoothly with increasing 
magnetic field, eventually saturating at high fields. Sim- 
ilar to the H\\c case, Lc shows the strongest effect, and 
La moves in the opposite direction. The magnetization 
for -ff||a is shown for comparison in Fig. [5] 

The critical fields Hci and Hc2 can be extracted from 
the magnetostriction data ALc for -ff||c. We take the 
peak in the second derivative of ALc{H) to be the anti- 
ferromagnetic transition. We have used this technique 
to map the phase diagram as a function of tempera- 
ture and field angle. The resulting temperature-field 
phase diagram is shown in Fig. [51 along with that deter- 
mined from specific heat and magnetocaloric effect data 
for comparison]^ The phase diagram from magnetostric- 
tion closely tracks that determined from previous mea- 
surements, however, the phase transitions are shifted to 
slightly lower fields. The discrepancy between the phase 
diagrams cannot be accounted for by errors in the mag- 
netic field centering, or by rotational misalignments. We 
assume therefore that this shift is due to pressure caused 
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FIG. 7: Critical fields Hci and Hc2 at T = 25 mK as a func- 
tion of field angle from the tetragonal c-axis determined from 
magnetostriction data. 



by the spring-mounted dilatometer. 

The critical fields Hd and Hc2 were also determined 
for different field angles 9 to the tetragonal c-axis. The 
field-angle phase diagram is shown in Fig. [7l The region 
of antiferromagnetic order between Hd and Hc2 shrinks 
as increases, with antiferromagnetic order finally van- 
ishing above 9 = 55 deg. 



III. CONCLUSIONS 

In conclusion, the compound NiCl2-4SC(NH2)2 shows 
significant magnetoelastic effects, which indicate a strong 
coupling between the magnetic spin system and the lat- 
tice. As discussed previously, the observed magnetostric- 
tion data can be accounted for by magnetic forces result- 
ing from the canted antiferromagnetic order, and by the 
Poisson effect that conserves the volume of the tetragonal 
crystal. The magnetostriction effects along the crystal- 
lographic c-axis are so large that the antiferromagnetic 
phase transition can be extracted from the data. The 
phase diagram as a function of temperature, field, and 
angle has been determined from the data, and agrees 



roughly with specific heat and magnctocaloric effect mea- 
surements. 

The original motivation for these measurements was 
to account for a discrepancy between spin-wave theory 
prediction of i?c2 = 10.85 T and the observed value of 
Hc2 = 12.6 T<^ The theory assumes constant values of 
Ja, Jci and D, whereas a magnetostriction effect would 
cause Ja, Jc, and D vary with field. 

However, this large discrepancy in Hc2 has been mostly 
resolved by taking into account spin fluctuations. Never- 
theless, small discrepancies between the theory and the 
experimental phase diagram and the magnetization vs 
field data still exist These small discrepancies might 
still be accounted for by the structural distortions of the 
sample with field. The theoretical predictions for Hd and 
Hc2 are a 'best fit', obtained by varying the measured pa- 
rameters D, Ja and Jc within their error bars. Thus both 
Hd and Hc2 show small discrepancies between theory 
and experiment. It is possible that a fit that preferen- 
tially matches Hd would result in a larger discrepancy in 
Hc2, which is in fact accounted for by magnetostriction 
effects at high fields. 

An important question that arises from these results 
is whether the tetragonal symmetry of the crystal is pre- 
served or whether the magnetostriction effects induce a 
orthorhombic distortion. If the tetragonal symmetry of 
the crystal is compromised in high fields and the energy 
scale of these distortions is significant then the Bose- 
Einstein Condensation description of the magnetic phase 
transitions might no longer be valid. It would only be an 
accurate description at Hd where the structural distor- 
tions go to zero. In order to investigate the possibility 
of structural distortions, resonant ultrasound measure- 
ments are in progress. 
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